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Abstract 

We compute the relation between the pole quark mass and the minimally sub- 
tracted quark mass in the framework of QCD applying dimensional reduction as 
a regularization scheme. Special emphasis is put on the evanescent couplings and 
the renormalization of the e-scalar mass. As a by-product we obtain the three-loop 
on-shell renormalization constants and Zg*^ in dimensional regularization and 
thus provide the first independent check of the analytical results computed several 
years ago. 
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1 Introduction 



In quantum chromodynamics (QCD), like in any other renormalizeable quantum field 
theory, it is crucial to specify the precise meaning of the parameters appearing in the 
underlying Lagrangian — in particular when higher order quantum corrections are con- 
sidered. The canonical choice for the couphng constant of QCD, as, is the so-called 
modified minimal subtraction (MS) scheme [1] which has the advantage that the beta 
function, ruling the scale dependence of the coupling, is mass-independent. On the other 
hand, for a heavy quark besides the MS scheme also other definitions are important — 
first and foremost the pole mass. Whereas the former definition is appropriate for those 
processes where the relevant energy scales are much larger than the quark mass the pole 
mass is the relevant definition for threshold processes. Thus it is important to have precise 
conversion formulae at hand in order to convert one definition into the other. 

By far the most loop calculations performed within QCD are based on dimensional regu- 
larization (DREG) in order to handle the infinities which occur in intermediate steps. It is 
well known, however, that it is not convenient to apply DREG to super symmetric theories 
since it introduces a mismatch between the numbers of fermionic and bosonic degrees of 
freedom in super-multiplets. In order to circumvent this problem and, at the same time, 
take over as many advantages as possible from DREG the regularization scheme dimen- 
sional reduction (DRED) has been invented (see, e.g., Refs. [2,3]). Indeed, the application 
of DRED to supersymmctric theories leads to a relatively small price one has to pay at the 
technical level. An elegant way is to introduce an additional scalar particle (the so-called 
e scalar) at the level of the Lagrangian and to proceed for the practical calculation of 
the Feynman diagrams as in DREG. The situation becomes more complicated in case the 
symmetry between fermions and bosons in the underlying theory is distorted, e.g., after 
some heavy squarks have been integrated out. In such situations couplings involving the e 
scalar, which are called evanescent couplings, renormalize differently from the gauge cou- 
plings and therefore one has to allow for new couplings in the theory. The same is true if 
DRED is applied to QCD: in addition to as four new couplings have to be introduced, each 
of which has their own beta function governing both the running and the renormalization. 
In this paper we take over the notation from [4] and denote them by ae (the coupling 
between e scalars and quarks) and rji {i = 1,2,3), which describe three different four-£ 
vertices. 

The one-loop relation between the MS and pole mass has been considered long ago in 
Ref . [5] . At the beginning of the nineties the mass relation to two-loop order [6] was one 
of the first applications of two-loop on-shell integrals. In a subsequent paper also the 
two- loop result for the on-shell wave function counterterm has been obtained [7]. The 
three-loop mass relation has been computed for the first time in Ref. [8, 9] where the off- 
shell fermion propagator has been considered for small and large external momenta. The 
on-shell quantities have been obtained with the help of a conformal mapping and Fade 
approximation. The numerical results of Ref. [8,9] have later been confirmed in Ref. [10] 
by an analytical on-shell calculation. The three-loop result for has been obtained in 
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Ref. [11]. In this paper we provide the first independent check of the analytical results 
for and in the DREG scheme. 

The renormalization scheme based on DRED together with modified minimal subtraction 
is called DR. As far as the relation between the DR and the pole mass is concerned one can 
find the one-loop result in Ref. [12]. The two-loop calculation [13] has been performed 
in DRED, identifying, however, the evanescent coupling with In this paper we will 
provide the more general result and furthermore extend the relation to three-loop order. 

The paper is organized as follows: In Section [2] we outline the general strategy, provide 
technical details and derive the result for the on-shell mass and wave function counterterm 
in the case of dimensional regularization. The peculiarities of dimensional reduction are 
explained in Section [3] where all relevant counterterms are listed. Finally, Section H] con- 
tains the main result of our paper: the relation between the pole mass and the minimally 
subtracted mass in the framework of dimensional reduction up to three-loop order. 



2 On-shell counterterms for mass and wave function 
in dimensional regularization 

In order to compute the counterterms for the quark mass and wave function one has to 
put certain requirements on the pole and the residual of the quark propagator. More 
precisely, in the on-shell scheme we demand that the quark two-point function has a zero 
at the position of the on-shell mass and that the residual of the propagator is —i. Thus, 
the renormalized quark propagator is given by 

■70s 

Spiq) = -, r (1) 



-i 



(2) 



where the renormalization constants are defined as 

m,,o = ^^'Mg, (3) 



^0 = ^JZ^^iJ. (4) 

ip is the quark field with mass rrig, Mg is the on-shell mass and bare quantities are denoted 
by a subscript 0. S denotes the quark self-energy contributions which can be decomposed 

as 

E(g,mg) = mgSi(g^mg) + -mg) S2(g^mg) . (5) 



The calculation outlined in Ref. [6] for the evaluation of and reduces all occurring 
Feynman diagrams to the evaluation of on-shell integrals at the bare mass scale. In 
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particular, it avoids the introduction of explicit counterterm diagrams. At three-loop order 
we find it more convenient to follow the more direct approach described in Ref. [10, 11], 
which requires the calculation of diagrams with mass counterterm insertion. Following 
the latter reference, we expand S around = 



_d 



d 



dq^ 



M,)\2M'—T.,{q\M,) ^ + ^^(M^^ M,) + . . . . (6) 



Inserting Eq. ([6]) into Eq. ([T]) and comparing to Eq. ([2]) we find the following formulae for 
the renormalization constants 



Z^' = 1 + Si(M2,M,), (7) 



g^=Aq - 



If we consider the external momentum of the quarks to be g = Q(l + t) with = M^, 
the self-energy can be written as 

E(g, M,) = M,Si(g2, M,) + {0 - M,)S2(g^ M,) + tQ E^iq^ M,) . (9) 
Let us now consider the quantity Trj^^^S} and expand to first order in t 

Tr|^^^S(g,M,)} = ^,{q\ M,) + t^,{q\ M,) 



d 



= Si(M2,M,) + (^2M2— Si(g2,M,) ^^^^^ +S2(M2, M,)J t 

+0 (t^) . (10) 

Thus, to obtain Z^^ one only needs to calculate Si for = M^. To calculate Z^^, one 
has to compute the first derivative of the self-energy diagrams. The mass renormalization 
is taken into account iteratively by calculating one- and two-loop diagrams with zero- 
momentum insertions. 

Sample diagrams contributing to the quark propagator are shown in Fig. [H All occurring 
Feynman integrals can be mapped onto J^^ and L^'* as given in Eq. (4) of Ref. [14] and 
to seven more similar ones which will be listed in Ref. [15]. 

All Feynman diagrams are generated with QGRAF [16] and the various topologies are iden- 
tified with the help of q2e and exp [17, 18]. In a next step the reduction of the various 
functions to so-called master integrals has to be achieved. For this step we use the so-called 
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id) (e) (/) 

Figure 1: Sample three-loop diagrams. Solid lines denote massive quarks with mass rriq 
and curly lines denote gluons. In the closed fermion loops all quark flavours have to be 
considered. 

^^-^ 

//\\ ^/"Z'x^ 

/ / \ \ " "^^ ^ - ^ ^ 
^-4 1 ^1-4 

(a) (b) (c) 

Figure 2: Three-loop master integrals. Solid lines denote massive and dashed lines are 
massless scalar propagators. 

Laporta method [19,20] which reduces the three-loop integrals to 19 master integrals. We 
use the implementation of Laporta's algorithm in the program Crusher [21]. It is written 
in C++ and uses GiNaC [22] for simple manipulations like taking derivatives of polynomial 
quantities. In the practical implementation of the Laporta algorithm one of the most 
time-consuming operations is the simplification of the coefficients appearing in front of 
the individual integrals. This task is performed with the help of Fermat [23] where a spe- 
cial interface has been used (see Ref. [24]). The main features of the implementation are 
the automated generation of the integration-by-parts (IBP) identities [25] and a complete 
symmetrization of the diagrams. 

The results for the master integrals can be found in Ref. [11]. We have checked the 
results numerically with the Mellin-Barnes method [26,27] and the program MB [28]. We 
did, however, find some differences to Ref. [11]: In addition to the 18 master integrals 
given in that reference, we found the integral depicted in Fig. M^a) where the result can 
be found in Eq. (A. 2) of Ref. [14]. Furthermore, we already pointed out in Ref. [14], that 
the result for Fig[2](6) (denoted by /ig in Ref. [11]) is wrong. As it turns out, the result 
given in [11] corresponds to the integral depicted in Fig. ^c). Since there is a one-to-one 
correspondence between the two integrals with regard to the IBP relations, it does not 
matter which is chosen to be the master integral. 
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Most of the master integrals of Ref. [11] were already calculated in Ref. [19]. However, it 
turns out that there are differences in the order O (e) terms of the integrals 12-/7 of these 
references. The expressions given in Ref. [11] agree with our numerical results. 

The results for the renormalization constants can be cast into the following form 

with i G {m, 2}. It is convenient to further decompose the three-loop contribution in 
terms of the different colour factors 

+CFTFni {Cp Zr' + Ca Z[^^ + T^m Z^^ + T^n, 

+CpT^n, {Cp Zr^ + Ca Z^" + TfUu Zf^^) , (12) 

where Cp — {N^ — l)/(2iVc) and Ca — Nc are the eigenvalues of the quadratic Casimir 
operators of the fundamental and adjoint representation of SU (Nc) , respectively. In the 
case of QCD we have A''^ = 3. Tp = 1/2 is the index of the fundamental representation 
and Uf = 111 + Uh is the number of quark flavours, rii and Uh arc the number of light and 
heavy quark flavours, respectively. The former are considered to be massless, while the 
latter have mass niq. Although we have = 1 in our applications, we keep a generic 
label which is useful when tracing the origin of the individual contributions. as{fj) is the 
strong coupling constant defined in the MS scheme with Uf active flavours. 

The one- and two-loop contributions to the mass renormalization constant are 
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where = ln(/i^/Mg). C„ is Riemann's zeta function with integer argument n and a. = 
Li4(l/2). We give the one- and two-loop results to order O {e^) and 0(e), respectively. 
In general these terms are necessary for three-loop calculations. 

The individual three-loop terms read 



7FFF 



9 



128e3 

9 

16^ 



63 27 , 

\ 

256 128 



189 



^16^^ 256^^ 256 ^ 
4 128^ 



480' 



1 



457 111 2 3 2, . 
512 512 8 



14225 6037 



-TT 



3072 3072 
1 2^ 5^ 



TT^ + 57r^ In 2 



-ln^2 



3a4 



+ 



1371 333 2 9 2,. 27, 

\ TT^ - -TT^ In 2 H Ca 

512 512 8 16^^ 



512 256 ^ 



(15) 



7 



33 / 49 33 \ 1 /3311 43 2 3 2, o 
"l28i^ + ^768 "l28^^y'^+ ^1536 "512'' +16'' ^""^ 

9 , 379 , 33 .\ 1 100247 6545 o 25 2, . 

Cs H ^« H L,, - H ^ TT^ H TT^ In 2 

32^ 256 ^ 256 V e 9216 9216 36 

89 2, 2^ 3995, 1867 4 19 2. 45, 35 , 35 

TT^ In^ 2 Ca + ttVs + —U In^ 2 04 

72 384 8640 16 16^^ 144 6 



+ 



/14311 1135 2 71 2, 71 \ , 1797 . , , 

(W - 1536^ + 2 - -C3 j L„ + + , (16) 



121 1679 11413 1322545 1955 , 115 n, 

+ 77^ ^ = 7-. - — TT^ In 2 



576e3 3456e2 20736e 124416 3456 72 

11 2, 2^ 1343, 179 4 51 2. 65, 11, . 11 

+ — TT^ In^ 2 + C3 + — 7r\3 Cs H In^ 2 + —04 

36 288 ^ 3456 64 32^^ 72 3 

/ 13243 11 2 11 2, ^ ll>\r 2341 „ 121 „ 

+ \- —T^ In 2 + — C3 -L„ Ll Ll , (17) 

V 1728 72 24 16^V 1152 " 576 ' ^ ' 



3 ^ \ 1 / 65 7 2 1^ 23^ 3 ^2^ 1 



32e3 ^\ 192 + 32^'^ J (,384 + 128"" + 4^' + 64^'' + 64 ^'^ 7 e 

575 1091 2 11 21 ^ 2 2, 2c. 145, 119 4 1 4^ 

H ^ TT^ TT^ In 2 + -TT^ In^ 2 + Cs + - In^ 2 

2304 2304 9 9 96 ^ 2160 9 

8 /497 5 2 1 21 ^ ^Ar 117,2 llrs 

+ -04 Ti + -TT^ In 2 + -Cs L„ L^, Lf, , (18) 

3 V384 384 3 4^V 128 ^ 64 ^ ' 



11 121 / 139 1,\ 1 70763 175 2 H 2, ^ 
+ + 7C3 - + ttttt; + -twt;'^ + ttt^^ 1^2 



72e3 432e2 ^1296 4^7 e 15552 432 18 
1 2, 2^ 89 , 19 4 1,4^ 4 

/869 7 2 1 2, ^ 1. \ . 373 . 11 . 



r 1 5 35 2353 13 2 7 , 



36e3 216e2 1296e 7776 108 18 



1 5 35 5917 13 2 2, 



18e3 108e2 648e 3888 108 9 
108 18 / 36 18 



"3 l.Ai Hi._ii3 pi) 



8 



7FFH 



7FAH 



+ 



32e3 
_5257 
~2304 " 
8 

+ 3a4- 
11 



192 



+ 1 - 

5 



281 17 2 1 > 23, 3 



384 128 



64 



91 



^^TT^ + — TT^ In 2 - -TT^ In^ 2 + — C3 H n' . 

6912 36 9 96^ 2160 9 



/1145 
121 



221 o 1 2, 1> 



117 n 

128 64 



+ 



/ 139 1, \ 1 
( + -C3 - 



72e3 432e2 ^1296 4 
1 2, 2^ 109, 



144959 449 2 32 „ 

TT^ + — TT^ In 2 

15552 144 9 

43 4 1 2. 5^ 1.4 



+ 



1080 

583 13 2 1 2, ^ ^Ar 373 ,2 

TT^ + -T^ In 2 + -Cs L„ + Ll + 

108 36 6 2^V '^ 288 



18 



(22) 



(23) 



7FHH 



1 



+ 



36e3 216e2 
197 1 . 

r 

216 9 



+ 



35 



1296e 
13 



9481 4 2 11. 
7776 ^135^ ^18^^ 



36 ^'^ 



(24) 



For the wave function renormalization constant, we have at the one-loop level 5^2^^ = 
5Zm ■ The two-loop contribution to the wave function renormalization constant reads 
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Again, we give the result to order O (e), which is necessary for three- loop calculations. 
The individual three-loop terms are given by 
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Starting from the three-loop level, the wave function renormalization constant depends 
on the gauge parameter, ^. The parameter in the above equations is defined through the 
gluon propagator as 

<.(fc) = -^(^.^-e^)^"> (36) 
where a and h are colour indices. 

We want to mention that our results for Z'^ and agree with the literature [10, 11]. 
Whereas in [10,11] they are expressed in terms of the bare coupling we decided to use the 
renormalized as an expansion parameter which to our opinion is more convenient in 
practical applications. 

The genuine three-loop integrals which appear in the DREG calculation are the same for 
DRED. The main complication is the more involved renormalization which is discussed in 
more detail in the next Section. 



3 Renormalization in DRED 

Let us in this Section collect the DRED counterterms needed for our calculation. In addition 
to the strong coupling also the evanescent coupling has to be renormalized to two- 
loop order. The evanescent coupling^ r/i, 772 and 773 appear for the first time at three- loop 
order and thus no renormalization is necessary. Both for the heavy quark mass, m^, 
and the e-scalar mass, m^, two-loop counterterms are necessary. Whereas the couplings 
are renormalized using minimal subtraction the masses are renormalized on-shell. The 
corresponding counterterms are defined through 

K)' = ^.^^S- (37) 

We attach an additional index "DR" to the quark mass renormalization constant in order 
to remind that it relates the pole mass to the bare mass in the DRED schem^ (in this 
context see also Ref. [30]). 

Recently the quantities Z^. and Z^^ have been computed to three- and four-loop or- 
der [4,30], respectively. The results have been presented in terms of the corresponding 
(3 functions. For completeness we present in the following the two-loop results for the 

^ We refer to [4, 29] for a precise definition of the evanescent couplings. 

^ For the SU(i) gauge group, which we exclusively consider in this paper, there are three independent 
such couplings [29]. 

•^In principle such an index would also be necessary in Section [2] However, since the MS scheme in 
connection with DREG constitutes the standard framework we refrain from introducing an additional index 
there. 
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Note that we set A^c = 3 in all terms containing the couplings rji since our implementation 
of these couplings is only valid for SU{3). We would also like to point out that the 
analytical form of is identical to the corresponding result in the MS scheme. This 
has been shown by an explicit calculation in Ref. [3]. The one-loop result for Zg can be 
found in Ref. [29]. 

The one-loop corrections to Z^^'^^ can be found in Ref. [12] and the two-loop terms have 
been computed in Ref. [13]. For our calculation we also need the ^(e^) and C(e) parts of 
the one- and two-loop terms, respectively. 

In Section H] we want to present the finite result obtained by considering the ratio 
Z^^'^^/Z^. The quantity Z^ has been computed in Ref. [30] to three and in Ref. [4] 
even to four-loop order. Whereas in [4, 30] only the anomalous dimensions are given we 
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want to present the explicit result for the renormalization constant 
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In order to achieve the finite result for the relation between the pole and the DR quark 
mass it is necessary to fix a renormalization scheme also for the mass of the e scalar, 
m^. Although, there is in general no tree-level term in the Lagrangian, there are loop 
induced contributions to which require the introduction of corresponding counterterms. 
The relevant Feynman diagrams contributing to the £-scalar propagator show quadratic 
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id) (e) if) 

Figure 3: One- and two- loop Feynman diagrams contributing to the e-scalar propagator. 
Dashed hues denote e scalars, curly lines denote gluons and solid lines denote massive 
quarks with mass rriq. 

divergences and therefore, one needs to consider only contributions from massive particles. 
Thus, in our case, only diagrams involving a massive quark have to be taken into account. 
Some sample diagrams are shown in Fig. [31 

It is common practice to renormalize on-shell and require that the renormalized mass is 
zero to each order in perturbation theory [31]. This scheme is known as the DR' scheme [31] 
and offers the advantage that the e-scalar mass completely decouples from the physical 
observables. For supersymmetric theories the DR and DR' renormalization schemes are the 
same, while for theories with broken supersymmetry the latter one is most convenient. At 
one- loop order there is only one relevant diagram (cf. Fig.[3t^a)) which has to be evaluated 
for vanishing external momentum. A closer look to the two-loop diagrams shows that 
they develop infra-red divergences in the limit (cf., e.g.. Fig. [3](e)). They can 

be regulated by introducing a small but non-vanishing mass for the e scalars. After the 
subsequent application of an asymptotic expansion [32] in the limit = -C 
the infra-red divergences manifest themselves as ln(m£) terms. Furthermore, one-loop 
diagrams like the ones in Fig. [3](6) and (c) do not vanish anymore and have to be taken 
into account as well. Although they are proportional to m^, after renormalization they 
induce two-loop contributions which are proportional to , partly multiplied by ln(m£) 
terms. It is interesting to note that in the sum of the genuine two-loop diagrams and 
the counterterm contributions the limit rrie — > can be taken which demonstrates the 
infra-red finiteness of the on-shell mass of the e scalar. 

Taking the infra-red finiteness for granted, it is also possible to choose = ml = from 
the very beginning. Then the individual diagrams are infra-red divergent, however, the 
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sum is not. We have performed the calculation both ways and checked that the final 
result is the same. It is given by 
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(41) 



where the constants are defined after Eq. ( |T4l) . The overall factor Uh in front of the one- 
and two-loop corrections shows that the renormalization of mg only infiuences those terms 
which contain a closed heavy quark loop. 



4 in dimensional reduction 

The approach to extract the on-shell mass counterterm in DRED can be taken over from 
DREG as described in Section [21 i.e. one considers the inverse quark propagator and 
requires that it has a zero at the position of the pole. Again, the counterterm diagrams 
are generated order-by-order in a generic way. 

The major complication as compared to the calculation in DREG is the appearance of 
the evanescent couplings and the e scalars. In particular, there are three different four- 
e vertices. This leads to many more Feynman diagrams which have to be considered. 
Whereas in the case of DREG about 130 diagrams contribute there are more than 1100 in 
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id) (e) (/) 

Figure 4: Sample three-loop diagrams contributing to the quark propagator which have 
to be considered additionally in case dimensional reduction is used for the regularization. 
Solid lines denote massive quarks with mass rrig, curly lines denote gluons and the e scalars 
are represented by dashed lines. In the closed fermion loops all quark flavours have to be 
considered. 



the case of DRED. Typical Feynman diagrams are shown in Fig. [T] and Fig. |H The more 
involved renormalization has already been discussed in Section [31 

We avoid to present the result for the divergent quantity Z^^'^^ but show the result for 
the ratio to the DR quantity which is finite. We cast our result in the following form 
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Since our implementation of the couplings r/i is only valid for SU{3) we furthermore set 
Nc = 3. Our one-, two- and three-loop results read: 
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For tte = C(s the two- loop result agrees with [13], while the three-loop one is new. 



There is a very strong cross check of the results in Eqs. (H3l) . (H4|) and (H5l) for the limit 
Uh = 0. The starting point is the relation between the MS and the on-shell mass in DREG 
which can easily be obtained from the results of Section [2] and the MS counterpart of 
Eq. (HOll (see, e.g., Ref. [33]). In this relation, which depends on a^^, both and 
are replaced by their DRED counterparts using Eqs. (4.2) and (4.3) of Ref. [4]. In this 
way we could verify the results for z^^'°^|n,,=o- This provides a strong consistency check 
both on the results presented in this paper but also on the approach used in [4] for the 
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extraction of the conversion formulae between the MS and DR quantities. 

In order to get an impression of the numerical size of the corrections, both in DREG and 
DRED, let us consider the relation between the minimally subtracted and the pole mass 

for the case of the bottom and top quark. As input we use af'^'^^{Mz) = 0.1189 [34], 
Mfc = 4.800 GeV and Mt = 171 A GeV [35]. 

In the framework of DREG it is straightforward to convert af'^'^^{Mz) to ai^'''^^(M{,) and 



a. 



(6), MS 



(Mt) using four-loop accuracj0 leading to 



bottom : mf^{Mb) = 3.953 GeV = 4.800(1 - 0.0929 - 0.0493 - 0.0342) GeV 

= (4.800 - 0.445 - 0.236 - 0. 164) GeV , (46) 

top : mp(Mi) = 161.1 GeV = 171.4(1 - 0.0461 - 0.0109 - 0.0033) GeV 

= (171.4-7.9-1.9-0.6) GeV, (47) 

with ai^^'^^{Mb) = 0.2188 and af^'^^(Mt) = 0.1085. The numbers given in the round 
brackets of the above equations indicate the contributions from the tree-level, one-, two- 
and three-loop conversion relation. 

Within DRED the numerical analysis gets more involved since four more couplings appear 
whose values are needed for /i = and /i = Mt. As mentioned in the Introduction, DRED 
is an appropriate scheme for supersymmetric theories where in the strong sector only one 
coupling constant is present — like in usual QCD using DREG. Since all supersymmetric 
particles are heavier than the electroweak scale it is necessary to match the full theory 
to the Standard Model at some properly chosen scale /i'^'^"^. At this step the additional 
couplings appear. 

For the computation of ai^\Mb) and ril^\Mb) we use /i"^*^^ = Mz, evaluate in a first step 
(5),DR^^^-j ^gjj-^g three-loop relation given in Ref. [4] and require 



a. 



'DR(M^) = a^^\Mz) = vPiMz) , 
4'\Mz) = V3\Mz) = . (4J 



In a second step the renormalization group functions, which are known to four {a 
three (oe) and one-loop order (rji) [4, 30] are used to obtain 

af^^{Mb) = 0.2241, 

af(M6) = 0.1721, 

rif\Mb) = 0.2152, 

r]!^\Mb) = -0.01798, 

r]i^\Mb) = -0.005777. (49) 



*We use RunDec [36] for the running and decoupling of as in the MS scheme. 
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In the case of the top quark we choose /i = Mj. It is necessary to know the couplings 
for six active flavours and thus we evaluate in a first step ai^^'^^(M^) = 0.1178 and pose 
the analogue requirements as in Eq. fHHl) with "(5)" replaced by "(6)" and Mz by M^. 
This leads to 

af'^{Mt) = af\Mt) = 7]f\Mt) = 0.1096 , 

ryf (M,)=r/f (M,) = 0. (50) 

Finally, we can insert the results of Eqs. ( H9l) and (|50l) into Eq. ( l42l) and obtain 

bottom : mp(Mfc) = 3.859 GeV = 4.800(1 - 0.1134 - 0.0495 - 0.033) GeV 

= (4.800 - 0.544 - 0.238 - 0.159) GeV , (51) 

top : mp(Mi) = 159.0 GeV = 171.4(1 - 0.0581 - 0.0105 - 0.0030) GeV 

= (171.4 - 10.0 - 1.8 - 0.5) GeV . (52) 

The perturbative expansion shows a similar behaviour as for the MS-on-shell mass rela- 
tion: the three-loop terms amount to 159 MeV and 500 MeV, respectively, and are thus 
far above the current uncertainty for the bottom quark mass [37] and much larger than 
the expected uncertainty for the top quark mass [38]. 



5 Conclusions 

The main result of this paper is the relation between the pole and DR quark mass to 
three-loop order in QCD where dimensional reduction has been used as a regularization 
scheme. The conversion formula has been obtained in analytical form where all occurring 
integrals have been reduced to a small set of master integrals with the help of the Laporta 
algorithm. Due to the occurrence of evanescent couplings when using DRED within QCD it 
is more advantageous to use in this case dimensional regularization. However, the latter 
cannot be used in supersymmetric models. Thus, our result constitutes an important 
preparation for similar calculations in supersymmetric extensions of the Standard Model. 

As a by-product we have obtained the corresponding relation and the on-shell wave func- 
tion renormalization using dimensional regularization. This constitutes the first check of 
the analytical results obtained about seven years ago. Furthermore, we can confirm the 
observation of Ref. [11] that depends on the gauge fixing parameter starting from 
three-loop order. 
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